Quantum spin Hall effect in a transition metal oxide Na2lr03 
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We study theoretically the electronic states in a 5d transition metal oxide Na2lr03, in which 
both the spin-orbit interaction and the electron correlation play crucial roles. Tight-binding model 
analysis together with the fisrt-principles band structure calculation predicts that this material is a 
layered quantum spin Hall system. Due to the electron correlation, an antiferromagnetic order first 
develops at the edge, and later inside the bulk at low temperatures. 



PACS numbers: 71.70.Ej, 75.30.Kz, 75.80.+q, 77.80.-e 



The non-trivial topology in condensed matter physics 
has attracted great interests over the decades as high- 
lighted by the celebrated discovery of the quantum 
Hall (QH) effect [l|, 0]. Haldane theoretically studied 
the QH effect on the honeycomb lattice even without 
the Landau levels 0], which suggested that the non- 
trivial topology is more ubiquitous in solids than ex- 
pected. A recent breakthrough in this field is the the- 
oretical and experimental discoveries of the quantum 
spin Hall (QSH) effect in time-reversal symmetric insu- 
lators 0, MM m i, i, [3, El, El, El m ■ intuitively it 
can be regarded as two copies of QH systems with up 
and down spins, but is driven by the spin-orbit interac- 
tion (SOI) instead of the external magnetic field. The 
Z2 topological number, which distinguishes a topologi- 
cal insulator, i.e. QSH insulator, from an ordinary band 
insulator [1] , is closely related to a Kramers doublet pro- 
tected by the time-reversal symmetry, and corresponds 
to the presence or absence of gapless helical edge modes 
in the semi- infinite system Q. The theoretical design of 
a topological insulator using HgTe/CdTe quantum wells 
[ist was followed by the experimental realization 14 1, 
and now this novel state of matter has been firmly estab- 
lished. 



However, topological insulators have been limited to 
semiconductors at low temperature. This is because we 
need the large SOI and the fine tuning of the band struc- 
ture. Therefore one important development is to realize 
more robust topological insulators at higher temperature 
by the larger SOI. Another interesting development is to 
study the interplay between the non-trivial topology and 
the electron correlation [l5| . Generally the electron cor- 
relation is stronger in d- and /-electrons than in s- and 
p-electrons. When we look at transition metal ions in the 
periodic table, the electron correlation is the strongest in 
3d elements and decreases to Ad and to 5d elements be- 
cause d-orbitals are more and more extended, while the 
SOI increases as the atomic number. Thus in 5d tran- 



sition metal oxides, both the SOI and the electron cor- 
relation become important with the same order of mag- 
nitudes. In addition, a variety of crystal structures and 
even tailor-made structures such as superlattice are avail- 
able in transition metal oxides. These advantages will be 
useful to design topological insulators. 

In this Letter, we study theoretically the electronic 
states of a newly synthesized compound Na2lr03 in 
terms of the tight-binding model analysis and the first- 
principles calculations as a representative example to pro- 
pose a way to design topological insulators in 5d transi- 
tion metal oxides by using the complex transfer integrals 
and the lattice geometry. This material is predicted to be 
(i) a QSH insulator, (ii) the edge antiferromagnet (AFM), 
and (iii) the bulk AFM with decreasing temperature. 

5d-orbitals are rather extended and subject to the 
large crystalline field. Under the octahedral crystalline 



field, d-orbitals are split into Cg {x — y , 3z 



'■)- 



and t2n {xy, yz, za;)-orbitals by lODq of the order of 
3eV |17| . The SOI is quenched in Cg-orbitals but remains 
effective in t2g-orbitals, which form effectively the triplet 
with ^cff = 1- Explicitly, {\yz) ±i\zx))/\/2 correspond to 
\£lff = ±1), while \xy) to \eig = 0). Including the SOI, 
we obtain the states with the total angular momentum 
joff = 3/2 and 1/2. The wavefunctions with jcs = 1/2 
read as 



I + 1/2) =(+|xj/ t) + \yz i) + i\zx 

I - 1/2) ={-\xy i) + \yz ]) - i\zx t))/V3. 



(1) 



The central idea is that the transfer integrals between 
these complex orbitals and oxygen orbitals become com- 
plex. For example, consider a p^-orbital. The transfer 
integral between | ± 1/2) and Pz is proportional to e^**, 
where 6 is the angle between the x axis and the bond 
direction. This complex transfer integral is responsible 
for topological states in iridates. Recently, a layered per- 
ovskite oxide Sr2lr04 was studied by angle-resolved pho- 
toemission (ARPES), X-ray absorption, optical conduc- 
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tivity, and first-principles calculations [13, . Ir^+ has 
five electrons, one of which is in a narrow band mainly 
composed by joff = 1/2-states described above, leading 
to a Mott insulator with the AFM order. These experi- 
ments confirmed that spin-orbit coupled jeS = 1/2-states 
are realized, even though Sr2lr04 itself is topologically 
trivial. 

Now we focus on Na2lr03 whose layered crystal struc- 
ture contains the honeycomb lattice as shown in Fig. [T] 
(a) (For the three dimensional structure, see Fig. 2(d).) 
Each Ir atom is surrounded by an octahedron of six O 
atoms, which leads to the similar energy level scheme as 
Sr2lr04, i.e. one electron in jcs = 1/2-states. There- 
fore we can construct the effective single-band model on 
the honeycomb lattice. Since the O p-level Cp are around 
3cV lower than the Ir d-lcvel [l3l , we can integrate out 
p-orbitals to obtain the following effective Hamiltonian 



Y,[4d, + h.c] + ^ \4i[^d, + h., 



where {ij) and {{ij)) denote the nearest-neighbor (NN) 
and next-ncarcst-ncighbor (NNN) pairs, respectively. 
The transfer integral t between a NN pair is real and 
spin-independent as given by 



3{pptt) 



3 ed 



(3) 



where (pdn), (ppcr) and (ppTr) are Slater-Koster parame- 
ters between pd and pp, respectively [2^. Note that the 
contributions of the order of {pdny /{ed — Cp) cancel out in 
the honeycomb lattice, in sharp contrast to Sr2lr04 with 
the square lattice. The transfer integral between a NNN 
pair depends on spin, leading to a topological insulator. 
The local x, y, and z axes at an Ir atom are chosen to 
point in the direction of neighboring O atoms as shown 
in Fig. [TJ Therefore Z ^ (.t + y + z)/%/3 is perpendicu- 
lar to the honeycomb plane. With this convention, the 
transfer integral is a 2 x 2 matrix in the spin space, and 
is written as 



(4) 



where a = x, y, z is the direction whose projection onto 
the honeycomb plane coincides with that of the hopping 
direction. The magnitude t' is given by 



6ed- ep 



[pdaf 



{pdaf 



V 



i^'d - <^p)(^d - Cp) 



(5) 



with e'j^ ^ ed + lODq. Note that the key to these complex 
transfer integrals is the asymmetry between two paths 
connecting a NNN pair. If there were an additional Ir 
atom in the center of the hexagon, leading to the trian- 
gular lattice, the transfer integral t' would vanish. The 
real transfer integral can be produced by the direct dd 




FIG. 1: (Color online) (a) The honeycomb lattice of Ir atoms 
in Na2lr03 viewed from the c axis. A large black circle shows 
an Ir atom surrounded by six O atoms (red small circles), (b) 
The transfer integrals on the honeycomb lattice. A black line 
shows —t, while blue dotted, red broken, and green dashed 
arrows indicate it'ax, it'o-y, it'az, respectively. 



hopping and breaks the particle-hole symmetry. How- 
(2) ever we put = for the moment since such term does 



not change the topological properties of the Bloch wave 
functions. 

To summarize these results, the transfer integrals are 
real and spin-independent for a NN pair, while complex 
and spin-dependent for a NNN pair. We can sec that 
this model is related to the Haldane model for the QH 
effect 3|^nd also to the Kane-Mele model for the QSH 
effect J,l5[. In fact, we can derive the effective Hamilto- 
nian around K and K' points of the Brillouin zone as 



(6) 

where ^{r) is the eight-component spinor field operator, 
and kx and ky are measured from K or K' points, ry's, 
r's, and cr's are the Pauli matrices for the valley {K or 
A"'), sublatticc (1 or 2), and spin (-1- or — ) degrees of 
freedom, respectively. This effective Hamiltonian is ac- 
companied with the spin Chern number 2 and belongs 
to the Z2 non-trivial class. However, two spin compo- 
nents are not decoupled away from K or K' points, and 
we need to analyze the lattice Hamiltonian Eq. ([2]) more 
seriously. With both time-reversal and inversion sym- 
metries, the Hamiltonian can be written in terms of five 
4x4 matrices ^a_= {tx, Ty, t^Gx, r^cTy, t^ct^) defined 
by Fu and Kane [9| as 



H^{k)=Y,da{k)Ta 



(7) 



a=l 



with di = — <(1 + cosfci -I- cosfc2), d2 = — 1(— sinfci + 
sinfc2), d^ = 2t'sinfci, di — 2t'sinfc2, and = 
—2t' sm{ki + ^2). We can apply the criterion by Fu 
and Kane and obtain the non-trivial Z2 number 
(—1)"^ = —1. In addition, we can directly find a pair 
of gaplcss helical edge modes by numerical diagonaliza- 
tion for the system with the open boundary condition in 
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one direction. The crossing point is at fc = tt for the 
zigzag geometry, while at fc = for the armchair geome- 
try, where k is the wavenumber along the edge. Anyway 
such crossing is protected by the Kramers theorem, and 
can get gapped only if the time-reversal symmetry is hvo- 
ken [lH. 

To justify the analysis above based on the tight- 
binding model and to estimate the transfer integrals 
quantitatively, we perform the first-principles band 
structure calculations by using the full-potential lin- 
earized augumented plane-wave method implemented in 



WIEN2k [2lj. Since the accurate crystal parameters of 



Na2lr03 are not yet known, we use those of Na2Pt03 and 
then perform atomic position relaxation. The mufhn- 
tin radii Rmt are set to 2.05, 2.10, and 1.65a.u. for Ir, 
Na, and O atoms, respectively. Wavefunctions inside 
atomic spheres are expanded in spherical harmomics up 
to ^max = 10, while those outside spheres are expanded in 
plane waves with the cut-off i?MT^max = 7.0. Figure [2] 
(a) shows the relativistic density of states (DOS) includ- 
ing the SOI, and (b) is its enlargement near the Fermi 
energy Ep = 0. First, the major part of the oxygen 2p- 
band is well below Ep ranging from —7 to — 3eV, but 
are hybridized with the 5d-band. The crystal field split- 
ting between Cg- and f2£(-orbitals is of the order of 4eV, 
which is much larger than that of 3d-orbitals. Therefore 
it is seen that the states near Ep is mainly composed of 
the i2g-orbitals hybridized with the oxygen 2p-orbitals, 
which supports the assumption of tight-binding analysis 
described above. 

As shown in Fig. [5] (b), these t2g-orbitals are further 
split by the SOI into joff = 1/2- and jeff = 3/2-states. 
From the center of gravity of each of these two partial 
DOS, the splitting is estimated of the order of leV, and 
the states near Ep is mostly those of jeff = 1/2-states. 
There is a finite DOS at Ep even though a dip is ob- 
served there. As shown in Fig. [2] (c), this is due to the 
overlap of the band dispersions of the conduction and 
valence bands. However, it is well known that the local 
density approximation (LDA) usually underestimates the 
band gap and it is plausible that the gap should open, 
considering the fact that the real material Na2lr03 is an 
insulator at room temperature above the magnetic phase 
transition at T/v = lOK. In Fig. [5] (c) shown the tight- 
binding fit (thick line) of the energy dispersions obtained 
by the first-principles band calculation (thin line). The 
tight-binding model includes the spin-independent real 
transfer integral tg between the NNN sites and the in- 
terlayer transfer t± besides the t and t' terms in Eq. ^ 
as shown in Fig. [2] (d). We set the transfer integrals 
to t = 310K, t' = lOOK, t'g = -130K, and t± = 60K, 
and the band structure obtained by the extended tight- 
binding model can roughly reproduce the features of the 
first-principles calculations. The deviations could be due 
to the mixing of Cg- and 2p-orbitals. For the present pur- 
pose, i.e., to study the topological nature of the Bloch 



wavefunction, this deviation does not matter. Especially, 
tg just shifts the energy dispersion and does not change 
the Bloch wavefunction itself. The helical edge chan- 
nels for a layer will turn to the two Dirac fermions at 
/ca = and tt on the surface of of the sample due to the 
non-zero interlayer coupling t±. Therefore this system 
belongs to a "weak" topological insulator in the classi- 
fication by Ref. Although the backward scattering 
between two Dirac fermions could induce localization, the 
"side" surface of the sample remains metallic as long as 
disorder is weak because this 2D system belongs to the 
symplctic class. Even in the localized case, the finite DOS 
on the surface including the c-axis which can be detected 
by scanning tunneling spectroscopy (STS), is a distinct 
feature from the usual insulator. Furthermore, the pres- 
ence of dislocation can hi ghli ght the non-triviality of a 
weak topological insulator [23|. To summarize, the tight- 
binding model analysis above is supported by the first- 
principles calculation provided that the band gap is larger 
than that of LDA estimation. 

Up to now we focused on the single-particle band 
structure, but in the real material the electron correla- 
tion energy U is of the order of 0.5eV as estimated in 
Sr2lr04 [18[. There arc two pictures for the effect of the 
electron correlation in the present case where two atoms 
are in the unit cell. One is the weak correlation case, 
where the band insulator can be antiferromagnetically 
ordered when U is larger than the critical value. The 
other limit is the Mott insulator, where the U localizes 
an electron at each site independent of the lattice struc- 
ture, and only the spin degrees of freedom survive. It 
is noted that these two pictures cannot be clearly dis- 
tinguished and cross over as U increases. Therefore it is 
important to study the two limits theoretically to under- 
stand the real material, which is most probably sitting in 
the intermediate region. 

In this Letter, we take the former view that the system 
is a topological band insulator, and the magnetic order- 
ing occurs due to the electron correlation at low temper- 
ature. As discussed above, the nontrivial topology leads 
to the helical edge modes at the boundary of the sample. 
Consequently, it is expected that the magnetic ordering 
occurs first along the boundary of the sample as the tem- 
perature is lowered, since the gapless helical edge chan- 
nels are more susceptible to the correlation. Because of 
the strong SOI, the spin anisotropy is expected. Within 
the mean field approximation, we have confirmed that the 
bulk AFM appears below Tn, and the edge AFM appears 
just above Tn, and both of these orders are within the 
Ary-plane, i.e., the easy-plane anisotropy. Considering 
the dispersion along the c-axis, the edge AFM moments 
are coupled between layers, leading to the AF-like spin 
model on the surface. Therefore we expect the Kosterlitz- 
Thouless transition above the bulk Tn- If a single layer of 
honeycomb lattice can be fabricated, we expect the novel 
one-dimensional correlated electrons along the edge. This 
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FIG. 2: (Color online) (a) and (b) The relativistic DOS in- 
cluding the SOI in two different ranges of energy. Black 
thick, red thin, green dashed, and blue dotted lines indicate Ir 
Jeff = 1/2-, Ir Jeff = 3/2-, Ir eg-, and O p-bands, respectively. 
The Fermi energy is set to zero, (c) The first-principles band 
structure (thin line) and the extended tight-binding model 
with typical parameters t = 310K, t' = lOOK, t'o = -130K, 
and t± = 60K (thick line), (d) The interlayer coupling t± is 
indicated by black dashed lines, while the other transfer inte- 
grals are shown in Fig. [T] (b) . Due to the monoclinic crystal 
structure, layers are not stacked in the simple way such as 
AB-stacked graphene. 



might be realized by a "plateau" on the top surface of 
the sample. Since the plateau edge is spatially separated 
from the edges of underlying layers, the interlayer cou- 
pling between the edges is reduced and gapless helical 
edge modes are robust. From the Tomonaga-Luttinger 
(XL) model analysis of the helical edge modes Uj, ll2], 
it is concluded that the electron-electron interaction pro- 
ducing the AFM order and the gap in the charge excita- 
tion is the Umklapp scattering accompanied by the spin 
flip processes. For this interaction to be relevant, the 
exponent K should be less than 1/2, which means the 



lie. Therefore the AFM ordering and the metal/insulator 
transition are always linked for the helical edge modes in 
the weak coupling analysis. Another interesting aspect 
is that the domain walls of the edge AFM have the frac- 
tional charge ±e/2 for each layer [2J|. It is expected that 
a minigap opening at the edge due to the AFM order can 
be detected by STS as temperature is lowered to ~ Tn, 
In conclusion, wc showed that the spin-orbit coupled 
d-orbitals in transition metal oxides lead to the topo- 
logically non-trivial electronic states. As an example, 
the newly synthesized compound Na2lr03 is found to be 
a QSH insulator in the paramagnetic phase within the 
tight-binding model as confirmed by first-principles cal- 
culations. Three-dimensionality, i.e. the interlayer cou- 
pling leads to a layered QSH insulator with two Dirac 
fermions. At low temperature, the electron correlation 
drives the system first to the surface AFM phase, then 
to the bulk AFM phase. Na2lr03 provides the first ex- 
perimental stage to study the effect of the electron cor- 
relation in QSH insulators. 
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Note added: After submission of this paper, we found 
a paper by Jackeli et al. which suggests that the spin ex- 
change interaction in the honeycomb lattice of Ir atoms 
leads to the Kitaev model [25[. Their approach from 
the strong correlation limit is complementary to our ap- 
proach from the weak correlation limit. 



strong repulsive interaction [ll|, [13]. For K > 1/2, there 



is no gap opening, and the edge channels remain metal- 
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